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Abstract 

Using the coordinate transformation method, we solve the one-dimensional Schrodinger equation 
with position-dependent mass(PDM). The explicit expressions for the potentials, energy eigenvalues 
and eigenfunctions of the systems are given. The eigenfunctions can be expressed in terms of the 
Jacobi, Hemite and generalized Laguerre polynomials. All potentials for these solvable systems 
have an extra term Vm which produced from the dependence of mass on the coordinate, compared 
with that for the systems of constant mass. The properties of Vm for several mass functions are 
discussed. 
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I. INTRODUCTION 



In recent years, the study of quantum systems with position-dependent mass 



become one of the active subjects 
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physics, semiconductor theory and other related fields. In the theoretical researches, many 
methods developed in the study of systems with constant mass have been generalized to 
the systems with PDM and a number of interesting results have been produced. Of those 
researches, the searching of exact solutions of the quantum wave equations (Schrodinger 
equation, Klein-Gordan equation and Dirac equation) for these systems is an important 
aspect. The factorization methodj^j, operator methods j^l, coordinate transformation 
25, Q, 

supersymmetric quantum mechanics 

Q Lie Algebra 
31 1, path integral approach 3^ and so on that exactly solve the quantum 
wave equations with constant mass have been extended to the systems with PDM. 

In the present work, we will use the coordinate transformation method to the PDM 
Schrodinger equations and obtain their solutions for several potentials. There are two 
kinds of coordinate transformations. The first one connects two different solvable poten- 
tials such as the Coulomb potential and harmonic oscillator, Rosen-Morse and generlized 
Scarf potentials [26]. The second one transforms the Schrodinger equation into the second 
order differential equation which has solutions of the special functions such as hypergeo- 
metric functions and confluent hypergeometric functions, and so may provide systematical 
studies to the exact solutions of the Schrodinger equation. When combined with the shape 
invariance, the coordinate transformation method provides an important way to classify the 
solvable potentials {25, 



28|. We will discuss the coordinate transformation method in the 



26 



second sense in this paper. 

The paper is organized as follows. In section |nj the coordinate transformation will be 
introduced and exact solutions of the one-dimensional Schrodinger equation with PDM for 
several potentials will be given. The eigenfunctions of these systems can be expressed in 
terms of Jacobi polynomials, Hemite polynomials and generalized Laguerre polynomials. In 
section IIIH the effective potentials produced from the dependence of mass on the position 
of the particle are given and their properties are discussed for several mass functions. In the 
last section, some remarks and discussions will be made. 
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II. COORDINATE TRANSFORMATION, PDM SCHRODINGER EQUATION 
AND THEIR SOLUTIONS 



When the mass of a particle depends on its position, the mass and momentum operator 
no longer commute, so there are several ways to define the kinetic energy operator of the 
quantum system |l| . In this paper, we will adopt the form of the kinetic energy introduced 
by BenDaniel and Duke'2], and the Hamiltonian with the position-dependent effective mass 
M(r) and potential energy V^(r) reads 



H 



+ V{r) = - 



V- 



m(r) 



V 



(1) 



2M(r) ■ ' 2mo 
where mo is a constant mass and m(r) is a dimensionless position-dependent mass. Us- 
ing natural units (mo = % = 1) and only considering one-dimensional system, we get the 
following Hamiltonian from Eq. (0) 



H 



d 



d 



dx mix) dx 



T^/ \ 1 rf^ m! d , 



2m dx"^ ' 2m'^ dx ' ^ " ^ 
where m' denotes the first derivative of m{x) with respect to the coordinate x. The PDM 
Schrodinger equation corresponding to the Hamiltonian Q reads 

d'^ip{x) m' dip{x 



+ 2m [E -V{x)]^/j{x) = 0. 



(3) 



dx'^ m dx 

Now, we make the following coordinate transformation to the eigenfunction iIj{x) in Eq. Q 
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^(x) = f{x)Figix)), (4) 

where F{g) is some special function and satisfies the following second order differential 
equation 

(5) 



Q{g) and R{g) in Eq. ^ are determined by the types of special functions. In this paper, 
we will consider three types of special functions. Substituting Eq. (jH) into Eq. Q, one gets 
the equation 



d^F 



+ 



9" 



+ 



2f' 



m 



dg^ [{g'y fa' mg'\ dg 



dF 



+ 



f" 



m'f 2m [E - V{x)] 



f{g'Y mf{g') 



+ 



(9') 



On comparing Eq. ^ with Eq. (jH)), we may have the following relations 

X g" 2/' m' 

[g'r fg ^g 



0. (6) 



(7) 



^ Jl_ _ ^'f , 2m [E - Vjx)] 



From Eqs. (|7j) and (jH)), we have 
E - V{x 



' Q^?' - ^ + — , (9) 



2m 

, 2 



R(g) ^ - -Q' 



+ ^ - G(m)] , (10) 



where G{z) = 7- — | (7-) • making integration to Eq. Q, we get the transformation 
function f{x) to be of the form 



/(x) = J-exp -/ Q{g)dg\, (11) 




which has an extra factor ^/m compared with that for the constant mass|25l|. 

Now, it follows that the PDM Schrodinger equation can be exactly solved if the forms of 
functions Q{g) and R{g) (or the types of the special functions F{g)) are given for a mass 
m(x) that is a function of space coordinate x, and if the right hand side of Eq. ()10|) can 
be divided into two parts: one depends on the quantum number n, but is independent of 
x, while another part depends on x rather than the quantum number n. For the above 
division, the former is the eigenvalue E of the system and the latter is the corresponding 
potential V{x). Such system is exactly solvable with eigenfunction ^/'(x), energy eigenvalue 
E and potential V{x). 

In the following sections, we will take F{g) to be of the Jacobi, Hermite and generalized 
Laguerre polynomials, respectively. And for the convenience, we introduce the auxiliary 
function 

/X 
^/m(x)dx. (12) 

G{g') in Eq.(|lU|) can be rewritten into the following form in terms of fi{x) 

Gig') = G(^') + (/^')'G(^), (13) 



where 

2 



^^dg_^ ^ f d9\ _ 3 

dfi dfi^ \dfi J 2 



d'^g^ f dg_ 



(14) 



which reduces to G{g') for the constant mass. Eqs. ()10|) and |T3|) show that there are a 
new term 4^(G'(/i') — G{m)) and an extra factor ^ in each term compared with that for the 



constant system |25l|. When m depends on the coordinate x, the term ■^{G{jj,') — G{m)) will 
be grouped into the potential V{x) of the system. In this sense, the PDM induces effective 
interaction between particles in the system, which is consistent with the concept of effective 
mass in the condensed matter physics and other related fields. We will denote Vm in the 
following sections as 
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4m 



1 



m 
m 



7 f m' 
4 V m 



(15) 



which is attributed to the dependence of the mass m on x. 



A. Jacobi polynomial and solvable potentials 



When we choose F{g) to be the Jacobi polynomial P^'^\g), the corresponding differential 



equation for P^'^\g) gives the expressions of Q{g)-, R{g) in Eq. (0) 
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Q{9) 



-a + (3 {2 + a + P)g 



1-9' 



1-9' 



(16a) 



Ri9) 



n{l + a + P + n) 



(16b) 



where a and f3 are parameters, n = 0, 1, 2 
1 



Substituting Eq. ^ into Eq. (HU)), we get 



E - V{x) 



9 



12 



,jGi9')-Gira)] + nil + a + P + n)-^^^^_^^^ 



+ [2{2 + a + (3) - {a - py] 



9 



12 



2\2 



99 



12 



^m(l — g 
+ [l-{l + a + pf] 



(17) 



9'9" 



4m(l -^2)2 ' L V ' ' J 8m(l -^2)2 

Now we will chose appropriate g{x) to make RHS of Eq. (fTTj) have a term that is independent 
of X , but may contain n. Once g{x) is determined, then we will get solvable potentials V{x) 
and its corresponding energy eigenvalue E from Eq. ()17|). Inserting g{x) into Eq. ()11|). we 
will obtain the explicit expression of f{x), and so the eigenfunction of the system is given 
by iIj{x) = f{x)Pn'^\g) for the above potential V{x). Similar to those for the systems with 
constant mass, there are two cases that satisfy the above requirements, and each has several 
different functions g{x). 

Case 1: g{x) satisfies the differential equation 



,/2 



g'^)m{x) 



C, 



(18) 



where C is a constant. In this case, one has 



-C 



1 



^9" 



(19) 



2(1-^/2) 

Each solution g[x) to Eq. (fTH|) corresponds to an exactly solvable system. According to 
the procedures given above, we can get its potential, eigenfunction and energy eigenvalues. 
We will not give the details of calculations and just list the results. Note that the energy 
eigenvalues are chosen such that -Eri=o = 0. The definitions of the parameters s, A are made 
so that the potential, eigenfunctions and energy eigenvalues can reduce to those for the 



system with constant mass when m{x) = 1, respectively 
(i) g{x) = i sinh(a/i(x)), C = —a? 



1 2 2 ^ 2t \2 

-s a a s — n , 



(20a) 



V{x) =-s'^a'^ + -a^(A^ — — s)sech^(a/i(a;)) 

a^A(2s + 1) tanh(a/i(x))sech(a/i(a;)) + Vm 



ip{x) = [m(x)]4 (cosh(ayu))~* exp |(— Atan~^[sinh(a/i(x))]} 

where s = — |(a + /? + 1), A = — a), 

(ii) g{x) = cosh(a/i(x)), C = —a? 

En = -s a - -a {s-n) , 



(20b) 



(20c) 



(21a) 



V{x) =-s'^a'^ + -a^(A^ + + s)cosech^(a/i(x)) 

— a^A ( s + - J coth(a/i(x))cosech(a/i(x)) + Vm 



1 — cosh(ayu(x)) ^ 2 
1 + cosh(ayu(x)) 



± 

ip{x) = [m(x)]* (sinh(ayu(x)))" 



where s = — |(a + P + 1), \ = — 1(/5 — a). 
(iii) g{x) = cos(a/i(x)), C = 

1 



1 



En = -^sV + -a2(s + n)2, 



(21b) 



(21c) 



(22a) 
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V{x) = — 2^'^^'^ + 2^'^^'^'^ + — s) csc^(a/i(x)) 
— -a^A(2s — 1) cot(a/i(x)) csc(a/i(x)) + 



(22b) 



1 — cos(a/i(x)) 

where s = |(a + /5 + 1), A = ^(/5 — a), 
(iv) g{x) = sin[a/i(x)], C = 



^(x) = [mix)]i (sinh(a/i(x)))^ ( ^ + cos(a/i(x)) \ ^ p^-^+-i^+-i)(cog(^^(3,))^ 



(22c) 



Er, 



1 2 2 , 1 2/ I \2 

■-S a H — a s + n) , 

2 2 ^ ' ' 



(23a) 



\^(x) = — 2'^^'^^ + 2*^^^"^^ + — s) sec^(a/i(x)) 
-Aa^(2s — 1) sec(a/i(x)) tan(a/i(x)) + 



(23b) 



^(x) = [m{xt^ (cos(a/x(x)))^ ( ! - sin(aMx)) ) ^ ^""'^'"'''^'"'^(sin(a/i(x))), (23c) 



where s = |(a + /5 + 1), A = |(/5 — a). 

Case 2: (^(a;) satisfies the differential equation 



(24) 



(1 — g'^ym{x) 
where C is a constant. Now, we have 

gA = -2C. (25) 

We will consider four solutions of g to Eq. (j23)- The corresponding potential, eigenfunction 
and energy eigenvalue of the exactly solvable system for each g are listed as follows: 
(i) g{x) = tanh(a/i(x)), C = 

A2 



1 X2 1 
771 -"-221^ 2 -"-2 

= -s a H -a a 

2 2s2 2 



(s — n)^ 



(26a) 



1 A^ 1 

\/(a;) = -s^a^ + —a^ - -a^s(s + l)sech^(a/i(x)) - Aa^ tanh(a/i(x)) + Vm, (26b) 

2 2^2 2 
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1 + taim[a(i{x)) J 

where 

a — s — n + a, P — s — n — a, a — . (27) 

s — n 

(ii) g{x) = coth(a/i(a;)), C — a? 

A2 



1 A2 1 

I— 1 f? 9 /\ Q J- O 

^" = 2" ^ + 2i^" - 2^ 



(s + ny + 



(s + n)^ 



(28a) 



1 1 

V{x) = -a^s"^ H -a? + -a^s(s - 1) csc^(a/x(x)) - Aa^ cot(a/x(x)) + Ki, (28b) 

a 

V'(a;) = [m{x)]^ (sinh(a//(a;)))"+^ ^coth(a^(a;)) - 1 V 



^coth(a/i(a;)) + 1/ (28c) 
X Pi-"-'^-"'-"-'*+^)(coth(a/^(x))), 

where 

a — —s — n-\-a, (3 — —s — n — a, a — . (29) 

s + n 

(hi) g{x) = —icot{an{x)), C = —a? 

= -V'^ + + - - (^' ^2°^) 

1 A^ 1 

\/(x) = — a^s^ H + -a^s{s + 1) csc^(a/i(a;)) - \c? coi{a\i{x)) + Ki, (30b) 

V'(x) = \m{x)f^ (sin(aAt(x)))"-"exp[aaAt(x)]P^^-"+^"'"-'*-^")(-icot(aAt(x))), (30c) 

where 

a = s — n^ia^ P = s — n — ia, a = . (31) 

s — n 

(iv) g{x) — — itan(a//(x)), C = — 

1 A^ 1 1 A^ 

= — a^s^ H -a^ + -0^(5 - nf - -a^- r^r, (32a) 

2 2s2 2 ^ ^ 2 (s-n)2' ^ ^ 

1 A^ 1 

^(3^) = ~ 2 "^^^ ^ 2^2 + 2 "^^^^ ^ sec^(a/i(x)) - Aa^ tan(a/i(a;)) + K., (32b) 
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^(x) = [m(x)]3 (cos(a/x(x)))""'exp[-aa/x(x)]Pi^-"+'^'^-"-"H-2tan(a/x(x))), (32c) 
where 

a = s — n + ia, P = s — n — ia, a = . (33) 

s — n 

It is obvious that all above results reduce to those for the systems with constant mass 
when m{x) = 1 and fi{x) = x From above explicit expressions for the potentials, 

eigenfunctions and energy eigenvalues of the solvable systems, we see that the energy eigen- 
values are the same as those for the systems with constant mass, but eigenfunctions and 
potentials do not so when the mass of the particle depends on x. The effects of PDM to 
eigenfunctions and potentials are twofold: the argument of function /i(x) and an extra factor 
containing m{x) in the eigenfunctions or a new term Vm in the potentials. In Section ITTH we 
will discuss the properties of Vm for several mass functions m{x). The above facts show that 
PDM will make the classes of solvable potentials more general than those for the constant 
mass. 



B. Hermite polynomial and solvable potentials 

When F{g) in Eq. (0) is the Hermite polynomial, i.e. F{g) = Hn{g), then Q{g) and R{g) 
in Eq. ^ have the following forms [33! Is^. respectively 

Q{g) = -2g, R{g) = 2n. (n = 0, 1, 2, ■ ■ ■ ) (34) 

With above Q{g), R{g) and Eq. (fTUI) . one has 

E - V{x) = ^ (2- + 1 - ^?^) + ^ m) - G{m)] . (35) 

There are two cases of conditions with which we can obtain the potentials and energy 
eigenvalues of the solvable systems from Eq. 
Case 1: If g{x) satisfies the equation 



m 

with > 0, we chose 



00, (36) 



g{x) = y/ujij,{x). (37) 
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Inserting Eq. (|H7|) into Eq. we get the eigenvalue and potential of the system, respec- 
tively 

En = nu!, (38a) 

V{X) = + [^(^)]2 ^ ^3g^^ 

From Eqs. (jS^), ()37|) and the transformation function reads 

/(x) = [m(x)]^exp(-^(72)^ (39) 



so the corresponding eigenfunction of the system with potential ()38b|) is 

i)n{x) = [m{x)]hxp{-^g^)Hn{g). (40) 

It is obvious that ()38a|) is the same as that for the harmonic oscillator up to a constant term. 
When m{x) = 1, Eqs. ()38b|) and pn)) reduce to the potential function and eigenfunction of 
the harmonic oscillator, respectively. 

Case 2: If g{x) is the solution of the differential equation 

where a; > 0, we have 

g{x) = ^f^^Hx)Y . (42) 

Putting Eq. (gSl) into Eq. dHHl), we get 

1/(1) = - I Ml)]-' - I + V„. (43b) 

Substituting Eq. ()42|) into Eq. one has the transformation function of the form 

/(x) = [m(x))u(a;)]^exp(-^/), (44) 

and the corresponding eigenfunction of the system reads 

i)n{x) = [m{x)fi{x)]Kxp{~g^)Hnig). (45) 
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C. Generalized Laguerre polynomial and solvable potentials 



When F{g) in Eq. (jS)) is the generahzed Laguerre polynomial, i.e. F{g) = L"{g), Q{g), 
R{g) in Eq. (0) will take the forms [ssl. 3^ . respectively 



Q{g) = 1, R{g) = -. (n = 0, 1, 2, ■ ■ ■ , a ^ -1, -2, -3, 

9 9 

Substituting Eq. into Eq. (fTTHl . we obtain the relation 



(46) 



Amg 



2mg^ 



-J\2 



1 



8m 4m 

Case 1: When g{x) satisfies the equation 

i9'y 



— [Gig') - G(m)] . 



4cj, 



mg 



with u; > 0, we have 



g{x) = u [fi{x)y 



(47) 



(48) 



(49) 



Inserting Eq. (jl^ into Eq. (jTTj) . we get the energy eigenvalues and potential of the system 



En = 2nu, 



(50a) 



V{x) 



l + ^]uj+}-uj^ [/i(x)] 



l2 , l/(/ + l) 



2 [^,ix)Y 



where 



I = a — 



2' 2' 2' 



In this case, the transformation function is 



(50b) 



(51) 



f{x) = [m{x)]-^ [n{x)]^+^ exp{-^ujn{xf). 
With Eqs. and ()52|1 . we obtain the corresponding eigenfunction of the system 
Mx) = [mix)]-^ [/i(x)]'+^exp(-^cu/i(x)2)L^^\a;/x(x)2). 



(52) 



(53) 



If / = 0, 1, 2, • • • , and / is viewed as the angular momentum quantum number, then Eq. ()5()a|l 
is the energy eigenvalues for the three-dimensional harmonic oscillator j35|. When m(x) = 1, 
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Eqs. fl50b|) and (fS^ reduce to the potential and eigenfunction for the three-dimensional 
isotropic harmonic oscillator, respectively. 
Case 2: If g{x) satisfies the equation 



(9 



l\2 



a\ (54) 



with a 7^ 0, we choose 

g{x) = exp[—afi{x)]. (55) 
With Eqs. (|^. and (|T7|). we obtain the energy eigenvalues, potential of the system 

En = ^a^s^ - {s - nf , (56a) 



where 



1 ^2 I 

V{x) = -ci^s^ + — exp [— 2a/i(a;)] — -a^exp [— a/i(a;)] + (56b) 

2 8 4 



1 In 
s = n+-a. (s^0,±-,--- ,±-,---) (57) 



'?/'n(a;) = [m(x)]'' exp [{n — s)a/x(x)] exp 



Lf-2")(exp(-a/i(x))). (58) 



From Eqs. (^^1 . (j^ . (jlj) and (fTT|) . we have the eigenfunction of the system 

1 

— ( 
2 

It is follows that the expression ()56a|) is the energy eigenvalue for the Morse potential with 
zero angular momentum. When m{x) = 1 and with the appropriate parameters a, b and c, 
then ()56b|) and (|58|) are Morse potential and eigenfunction with zero angular momentum, 
respectively. 

Case 3: If g{x) is the solution to the differential equation 



(9 



i\2 



= Au\ (59) 

m 

where a; > 0, we take 

g{x) = 2ujfi{x). (60) 

If we make the replacements u = , a = 21 + 1 {I ^ —1, — |, — |, ■ ■ ■ ) and use Eqs. 
and (jUj), we have 

" " 2(/ + l)2 " 2(n + / + l)2' 
12 



In this case, the transformation function can be written as 

1 . 
/(x) = [m(a;)]4 {fi{x)y^ exp 

so the eigenfunction of the system reads as 

V'n(a^) = (/i(x))'"'""^ exp 

If a = Z{Z is the charge numbers of the particle), / = 0, 1, 2, ■ ■ ■ , and / is regarded as the 
angular momentum quantum number, then Eq. (jfjlaj) is just the energy eigenvalues for the 
three-dimensional Coulomb potentialj^^. When m(x) = 1, Eqs. ()61b|) and (jHSI) reduce to 
the three-dimensional Coulomb potential and its eigenfunction, respectively. 



n + l + 1 



(62) 



n + / + 1 



/i(x) 



2^(2/+!) ( 



2a 



(63) 



III. MASS FUNCTIONS AND POTENTIALS 

In this section, we discuss the effective potentials Vm and their properties due to the 
dependence of mass on the coordinate for several mass functions. 

Example 1: We take the effective mass function to be of the form 

which has been used in many studies 0, Substituting Eq. into Eq. (fTT|. we get 

/i(x) = X + {b — 1) arctanx. (65) 

With Eq. ()15p. the contribution to the potential from mass function is 

jb - l)[3x^ + 2(2 - b)x' - b] 

2{b + xY ■ ^ ^ 

It is seen that Vm = when b = 1, which corresponds to the system of constant mass. 
When 0<6<lorl<6<4, Vm has three extreme points 



x = 0, x = ± I 6- 1 + yi(262-26 + 3)j . (67) 

If 6 > 4, then there are five extreme points for Vm, three of which has the same form as that 
in Eq. (jU7j) . the other two extreme points are 



a; = ± ( 6-1- -\/i(262- 26 + 3) ) . (68) 
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The characteristic curves for Vm with different values of the parameter b are depicted in 
Fig.l. It is seen that Vm behaves hke a barrier and it decrease as the parameter approaches 
to 1 from 6 < 1. While when b > 1, Vm just looks like a well. In this sense, Vm will bound 
the motion of the particle. 




FIG. 1: The potential Vm for the mass function (|64j) and the different values of the parameter 
b: (a) b = 0.5(solid line), b = 0.6(dashed line), b = 0.7(dot-dashed line); (b) b = 2.5(solid line), 
b = 3.0(dashed line), b = 3.5(dot-dashed line); (c) b = 9.0(solid line), b = 10.0(dashed line), 
b = 11.0(dot-dashed line). 



Example 2: The effective mass isjOL llOl l20| 



mix] 



(69) 



where 6 > to assure that mass is finite when x — > ±oo. Now, inserting Eq. ()fi9j) into (jllj) . 
we get 



-e 2 
b 



fi{x) = < 



2 i.. 

-62' 


X. 



The extra interaction due to the mass is 



Vrr, 



32 



(6^0, x>0) 
(67^0, x<0) 
(6 = 0) 



(70) 



(71) 



which is monotonously changed as |x| increases from 0. This Vm has the typical character- 
ization of Fig. 2 for various values of the parameter b. It is a kind of barrier whose width 
decreases with increasing the parameter b. 
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FIG. 2: The potentials Vm for the mass function 1)69^ and the parameters h = 0.4(sohd hne), 
h = 0.7(dashed hne), b = 1.0(dot-dashed hne), respectively. 



Example 3: The effective mass is of the form 11 1 

. . 1 



m[x] 



where 6 > to avoid its singularity. Now, we have 



(72) 



The effective potential Vm is 



/i(x) = ln(x -|- vh + x^) 

* m. 



(73) 



(74) 



8(6 + x2)' 

This Vm has the behavior of Fig. 3 for various values of the parameter h. It is a typically 
potential well whose width increases as the parameter h is increased. 
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FIG. 3: The potential Vm for the mass function (|72() and the parameters b = 2 (solid line), h 
4(dashed line), h = 6(dot-dashed hne), respectively. 
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Example 4: The effective mass is as follows 

m 

m{x) = 1 + tanh(6x), 

where a is a real parameter. Inserting Eq. (j75|) into Eq. (jllj) . one has 

V2 



^ tanh Y ^ (1 + tanh(6x)), [b ^ 0) 



X. {h = 0) 

Similarly, we get the potential produced from the dependence of mass on x 



32 



6^sech(aa;) [7cosh(6x) + sinh(6x)] [cosh(26x) — sinh(26x)] . 



(75) 



(76) 



(77) 



The curves of Vm for various values of the parameter b are displayed in Fig. 4. It can be 
seen that these Vm look like the semi-infinite potential barriers whose widths decrease as the 
parameter b increases. 




-1 

-1 . 5 
-2 

-2 . 5 



FIG. 4: The potential Vm for the mass function (jTSf) and the parameters b = 0.15(solid line), 
b = 0.2(dashed line), b = 0.3(dot-dashed line), respectively. 

With all above mass functions, we see that the dependence of the mass on the position 
of the particle will affect the behaviors of the system through two ways: the argument of 
the part in the potential that has the same form as that for the system with constant mass 
is fi{x) instead of x and an extra term Vm- These effects also increase the number of the 
solvable potentials for the same energy eigenvalue. 



IV. REMARKS AND DISCUSSIONS 

In this paper, we use the coordinate transformation method to study the exact solutions 
of the PDM Schrodinger equation for several potentials. The eigenstates of all these systems 
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with PDM can be expressed in terms of three kinds of special functions. We also give the 
explicit expressions of the potentials Vm for several mass functions that are used in some 
physically interested problems and study their properties. All these results will reduce to 
those for the systems with constant mass (see, for example 0, 26, |^|) if we set the mass 
functions m(x) to be constants in all equations in the above sections. It should be noted 
that there are two functions m{x) and g{x) in the PDM case, while only one function g{x) 
is concerned in the cases of constant mass. So, the classes with the same form of the energy 
eigenvalues but with different potentials related by coordinate transformation are enlarged. 

In our above discussions, we use the Hamiltonian of the symmetric form (0). If we adopt 
the Hamiltonian yi 



(78) 



with the condition a + /5 + 7 = —1, then all the results above also hold provided that we 
replace V{x) by K// in both Eq. (P) and other related relations, here 

1 TTl " 7T? ' 

V,ff{x) = V{x) + - [«(« + P+l) + f3 + l]- 

m 



.'2 



(79) 



' — and m" = This is so due to the fact that the Hamiltonian (f7S|) can be 



dx 

rewritten as 



dx^ 



H 



d 



d 



eff\ 



(80) 



^dx m{x) dx 
for the one- dimensional system. 

Also, for the three-dimensional systems with PDM and spherical symmetry, the solution of 

the system can be written as the product of angular and radial parts. The radial Schrodinger 

equation for the Hamiltonian takes the form 



where m' = 



1 

dm{r) 



d 1 d 

dr m dr 



0(r) + 



V{r) 



1 /(/ + 1) 



m 1 
2ni? r 



(f){r) = E(j){r), 



(81) 



2m 

is the radial wave function, E and I are the energy eigenvalue 



dr'^m- . 

and angular momentum quantum number of the system, respectively. Eq. ()81|) has the same 
form with the Schrodinger equation for the Hamiltonian (jHUj) . so the results in this paper can 
also be applied to the spherically symmetrical systems with PDM upon some modifications. 

We known th_at the special functions can have some generalizations, such as the q- 
deformed forms 34], so the corresponding differential equations have more general forms 
than Eq. (0). The coordinate transformation method can in principle apply to this gener- 
alized case and may give a more general classification to the solvable potentials. 
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